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Abstract Modelling Quark-Gluon Plasma formation and decay in high energy heavy ion 
reactions is presented in a framework of a multi-module setup. The collective 
features, governing the equlibrated fluid dynamical stages of the model are em- 
phasized. Flow effects formed from the initial conditions are discussed. Particu- 
lar attention is given to the improvement of the final hadronization and freeze-out 
part of the reaction which has strong effects on the observables. 



1. Introduction 

Phase Transitions imply that we have different phases of the matter we deal 
with, and it has an Equation of State (EoS). We can have an EoS if and only 
if the matter is in statistical equilibrium, or at least close to it. When this 
condition is satisfied we can apply a fluid- and thermodynamical description 
for the dynamics of the system. In high energy heavy ion reactions we observe 
many thousand particles produced in the reaction, so, we have all reasons to 
assume that in a good part of the reaction the conditions of the local equilibrium 
and continuum like behaviour are satisfied. 

The initial and final stages are, on the other hand, obviously not in statistically 
equilibrated states, and must be described separately, in other theoretical ap- 
proaches. The different approaches, as they describe different space-time (ST) 
domains and the corresponding approaches can be matched to each other across 
ST hyper-surfaces or across some transitional layers or fronts. The choice of 
realistic models at each stage of the collision, as well as the correct coupling 



2 



of the different stages or calculational modules are vital for a reliable reaction 
model. The most basic requirements are that conserved quantities should stay 
conserved when modules are interfaced to each other, and entropy should not 
decrease. 

For the theoretical description and understanding of high energy heavy ion 
colUsions it is important to mention that while our system is sufficiently large 
to assume statistical equilibrium, it is not a macroscopic system with Avogadro 
number of particles. We deal with a small, mesoscopic system, so discontinu- 
ities of the thermodynamical quantities, will not be observable, but at the same 
time the changes are sudden and rapid enough that they can be recognized. 

The other basic feature of mesoscopic systems is also important to us: in such 
systems fluctuations are not negligible but dominant, these can be observed and 
these may help to identify the properties of the phase transition. 

Dividing the reaction model into standardized modules and interfaces leads 
to additional advantages. Modules can be replaced easily, research groups can 
collaborate better and the GRIDS of High Performance Computers can increase 
the speed of complete reaction simulation by about an order of magnitude. 

2. Initial Stage Module 

Many fluid dynamical approaches use simple parametrizations for the initial 
state of the hydrodynamic stage with a few parameters. These can then be fitted 
or adjusted by comparing experimental data and model predictions. 

A more physical approach is to start from a dynamical, pre-equilibrium 
model, e.g. Parton Cascade Model (PCM), or Collective (or Coherent) Yang- 
Mills model, or different versions of flux- tube models. 

Our recent works are based on a latter type of approach [1], which was im- 
plemented in a Fire-Streak geometry streak by streak and upgraded to satisfy 
energy, momentum and baryon charge conservations exactly at given finite ener- 
gies [2, 3]. The effective string-tension was different for each streak, stretching 
in the beam direction, so that central streaks with more color (and baryon) 
charges at their two ends had bigger string-tension and expanded less, than 
peripheral streaks. The expansion of the streaks was assumed to last until the 
expansion has stopped. Yo-yo motion, as known from the Lund-model was not 
assumed. 

In the 1st version of the model at full stopping of string expansion, the 
baryon charge was distributed uniformly along each streak. This created an 
initial state for the hydrodynamic stage where each streak was stopped in its 
own CM frame, and had uniform baryon and energy density distribution. The 
sharp forward and backward end led, however, to a development of a large, 
Fwd/Bwd density peaks. [4] 
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Thus, we concluded that the assumption of uniform distribution and sharp 
cuts at the ends is overly simplified and has to be corrected. So, a Fwd/Bwd 
expansion was added to the model [3], which smoothed the Fwd/Bwd ends, and 
ehminated the artificial large, Fwd/Bwd density peaks. 

Our calculations show that such a tilted initial state leads to the creation of 
of the third flow component [5], peaking at rapidities \y\ « 0.75 [4]. Recent, 
STAR vi data [6] indicate that our assumption that the string expansion lasts 
until full stopping of each streak, may also be too simple and local equilibration 
may be achieved earlier, i.e. before the full uniform stopping of a streak. We did 
not explicitly calculate dissipative processes, some friction within and among 
the expanding streaks is certainly present and experiments seem to indicate that 
this friction is stronger. 

3. Relativistic Fluid Dynamics 

When parts of our system reach local statistical equilibrium in the ST, fol- 
lowing a pre-equilibrium development, we can describe our system using Com- 
putational FLuid Dynamics (CFD) and an Equation of State. The equations of 
Fluid Dynamics (FD) can be most simply derived from the Boltzmann Transport 
Equation (BTE). 

The relativistic BTE below, describes the time evolution of the single particle 
distribution function based on the assumptions that just two particle collisions 
are considered (binary collisions), the number of binary collisions at x is pro- 
portional to f{x,pi)x f{x,p2) and f{x,p) is a smoothly varying function 
compared to mean free path. So, one can show that BTE has the form [see 
Chapter 3 of ref. [7]] : 



where C is the colhsion integral, and k and I stand for different particle 
species or particle components. Then, using microscopic conservation laws 
in the collision integral, from the BTE we can derive the differential form of 
conservation laws [see section 3.6 of ref. [7]]: 
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where denotes d/d/j, and there is a summation for indexes occurring twice. 
These equations are vaUd for any distribution satisfying the BTE. These equa- 
tions can be used if we know the solution and then we can evaluate the conserved 
quantities. 

3.1 Perfect and Viscous Fluids 

Equations (1-3) are also the equations of FD. But these equations must not 
be considered as consequences of BTE, rather these equations are postulated, 
as the differential forms of conservation laws. This can be done because such 
conservation laws can be derived from other theoretical approaches also. 

These equations are not a closed set, because the energy momentum tensor 
and the particle four current should also be defined. In the Eulerian or perfect 
fluid dynamics we postulate the form of energy-momentum tensor, and an EoS, 
P = P{e, n), must be given. This provides a closed set of partial differential 
equations that is solvable. Any EoS, which is consistent with thermodynamics, 
can be used. One does not have to assume dilute systems with binary collisions, 
or the assumption of molecular chaos. The only requirement is the existence 
of local equilibrium, otherwise we would not have an EoS. 

In the case of Navier-Stokes fluid dynamics, we assume that the energy 
momentum tensor also contains the dissipative part r^'^(^), which allows for 
small, first order deviations from local equiUbrium. Thus, to have a solution 
we need not only the EoS, but also the transport coefficients that occur in the 
dissipative part of the energy-momentum tensor. 

Boltzmann's H-theorem impUes that if irreversible processes are present the 
entropy increases. In equihbrium the distribution is constant, the entropy is 
also constant, but it reached its maximum value. 

It is important to know that this can be seen from fluid dynamics also. In 
perfect fluid dynamics the flow is adiabatic, no entropy is produced. This can 
be proven using the equations of perfect fluid dynamics and standard thermody- 
namical relations. Entropy production is strictly, and quantitatively connected 
to dissipative processes. Perfect flow is adiabatic even in the presence of phase 
transitions, if the two phases are in thermal, mechanical (P), chemical and phase 
equilibrium all the time! On the other hand if the phase transition deviates from 
equilibrium this leads to entropy production. See Assignment 9.4.a in ref. [7]. 

As shock waves, detonations, deflagrations, are ideaUzations of very sharp 
or very rapid processes as discontinuities, the above conclusions are valid for 
discontinuities also. Nevertheless, frequently, dissipative processes can be ne- 
glected for most of the flow except the sharp or rapid changes. Thus, as an end 
effect the dissipative processes are frequently localized in sharp fronts or in hy- 
persurfaces of discontinuities. One should not, however, forget that dissipation 
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is due to transport properties and characteristic constants of phase transition 
dynamics, even if it happens in "discontinuities". 

3.2 Equations of Perfect Fluid Dynamics 

We assume that our system is not homogeneous, but the gradients are small, 
so local distribution can be an equihbrium distribution, e.g. a Jiittner distri- 
bution. Now the local n(x), e(x), s{x) and P = P{n, e) are known, and we 
assume that in LR, T'*^ is diagonal. Now using the conservation laws with 
introducing the apparent density as: 

N = wy , 

the continuity equation takes the famiUar form 

{dt + V grad) N = -M div v , 
similarly introducing 

M = = w -I'^Y (4) 
e ^ r00 = (e + Pv2)72, (5) 

the energy and momentum conservation will take the form 

{dt + V grad)M = -M(divv) - grad P (6) 
{dt + V grad)£: = — e divv — div(Pv) . (7) 

Thus to solve the equations of relativistic fluid dynamics we have to solve 
the same partial differential equations as in the non-relativistic case. We cannot 
immediately apply the EoS to obtain the solution, as the EoS is given in terms of 
invariant scalars. So, we have to solve in addition the set of algebraic equations 
above, which coimect N, M, and e to n, e, P and v, for every fluid cell at every 
time step. 

Important to note that the relativistic treatment is necessary not only at high 
velocities (i.e. if v ~ c) but also at low velocities when the pressure is not neg- 
ligible compared to the energy density (i.e., P e)a.sw = e + P appears in the 
relativistic Euler equation! This is frequently the situation for ultra-relativistic 
gases, hke the Stefan-Boltzmann or photon gas (or radiation pressure) where 
P = e/3. This is actually the indication of the fact that the constituents of the 
matter move with velocity c or close to it. 

In principle the perfect fluid dynamics is absolutely unstable. Any sponta- 
neous deviation from an exact solution, will grow exponentially, as their energy 
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Figure 1. The pressure, in units of [GeV/fm'^], in the reaction plane at a late stage (13 fm/c) 
of a 65+65 Gev/nucl. Au+Au collision, at b=0.5 b-max, using the Bag Model EoS for QGP. The 
outer, supercooled regions have negative pressure, while in the very center it is still large. The 
directed flow, is visible as the outside, negative pressure regions in the upper-right and lower-left 
comers, close to the beam, Z—, direction. The third flow component [5] is shown by a smaller 
negative pressure peak at the lower side of the plot, slightly right of the middle. The section of 
the reaction plane, [x,z], shown is 12.1 x 9.6 fm. 



cannot be dissipated away [8]. In other words the Reynolds number tends to 
infinity and the system tends to turbulence. Viscosity is needed to stabilize the 
flow. 

There are several numerical methods to obtain a solution. In all cases the 
solution is discretized in some way and a coarse graining is introduced. This 
coarse graining automatically leads to dissipative and transport properties, and 
the parameters of these transport processes can be determined by numerical 
experiments. The result is that all CFD models are dissipative or viscous, which 
is an advantage as the these are usually stable. By decreasing the cell size, the 
dissipative effects can be decreased in CFD, to a limit when instabilities, or 
turbulence occurs. As a side effect, CFD models lead to entropy production, 
as part of the kinetic energy is dissipated away due to the coarse graining or 
smoothing. 
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4. Kinetic Theory of Particle Formation, Escape and 
Freeze-out 

Most frequently relativistic kinetic theory describes the time, t, development 
of the single particle distribution function, f{p;x) = f{p,r,t), of conserved 
particles in the 6-dimensional Phase-Space (PS), [p, r] and time. Then the total 
particle number does not change with time. We assume (usually) that the par- 
ticles are on mass shell, so the 4-momentum has 3 independent variables. The 
phase space distribution is a density in the 6 dimensional PS and it is an invariant 
scalar, because a Lorentz boost contracts the system in the configuration space 
but elongates it in the momentum space. 

If we want to describe a particle source (or drain) from a matter, the total 
number of newly crated particles can be given as 



where g{p.. x) is an invariant scalar density in the 4-dimensional PS and 4- 
dimensional space-time (ST). If the particles are on mass shell the momentum 
space integral can be reduced to 



The particle formation or source density h{p, x) from a matter, can be char- 
acterized by 3 factors: 
PS distribution, f{p,r,t), 
PS emission probabihty, Pcsc-ps{p', x), 
ST emission probability, Pesc-ST{x;p). 

Here / can be an arbitrary distribution, but we usually discuss matter in (or 
close to) statistical equilibrium, with known PS distribution. Pesc-PS is most 
frequently assumed to be unity, while P^sc-ST is the ST distribution if the 
source (e.g. a 4-dimensionaI Gaussian in the ST). The experimental analysis 
of two particle correlations from heavy ion reactions, is based on these very 
simple assumptions in most cases. 

These are very strong assumptions. It is not supported by experiments that 
the emission probability is isotrope and homogeneous in the PS, or neither that 
the ST distribution of the source is Gaussian. 

4.1 Non-isotropic particle sources 

Not only in heavy ion reactions but in any dynamical process, particle cre- 
ation, (or condensation) happens mostly in a directed way: the phenomenon 





(8) 



where h{p, x) is also an invariant scalar density, as^^is also one. 
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propagates into some direction, i.e. it happens in some (space-like) layer or 
front (like detonations, deflagrations, shocks, condensation waves or freeze- 
out). The reason is that neighbouring regions in the front may interact, mini- 
mize the energy of the front by evening it out, provide energy to neigbouring 
regions to exceed the threshold conditions. Even in relativistic processes, that 
are time-like, (have time-Uke normal) and the neighbouring points of a front 
cannot be in causal connection, in the LR frame of the matter the dynamical 
process may and frequently has a direction. (See the example in ref. [9].) 

These fronts or layers are not necessarily narrow, but they have a character- 
istic direction (or normal, d^a^^), thus the ST integral, df'x, can be converted 
into 

d^x — > ds'' d^a^. 

The front can actually be directed time-like, d^a^^ d^a^ = +1, or space-like, 

dV^ dV^ = -1. 

An example for emission distributed in the ST, but being time-directed can 
be described by the expression: 

Nnew = j d^pj drd^'^^ feq{p,x) 9^e{p-d^a,) G {t - tfo{x)) (9) 



where 



„ eip'-d'^a,) (10) 

is the PS emission probabihty, where the step function eliminates the possibility 
of emission of particles, in the direction opposite to the normal of the front (for 
time-like normals it is always -1-1), and the term, ^ = (1, v), together with 
the 6- dimensional PS density, f^q, yield the generalized 4-dimensional flux of 
particles with momentum p: feqip, x) 

Here G is a Gaussian time distribution or emission density for time-like 
directed fronts: 



dr G{t — TFo) = — exp 



27r Tcoll. 



(r - TFof 
'^^coll. 



and for space-like fronts: 



ds G{s — spo) = ^ exp 

m.f.p. 



27r A, 



(s - SFO 



i2 



In more complex models the emission probabilities may take more realistic 
and more complicated forms. Furthermore, /e.,^. can be space-time dependent, 
and can be determined self-consistently during the detonation, deflagration or 
freeze out process [10, 11, 12, 13] . 



Phase Transitions in . 



9 



Let us rewrite the escape probability, (10), with the hyper-surface element, 
covering both the timelike and spaceUke parts of the freeze-out surface, 

n*.c-P5 = 1^ e^/^dV,). (11) 

Here and are ST dependent and this yields the secondary 

X— dependence of P*sc-PsiP'i ^^^^ velocity equal to Un = 

(1, 0, 0, 0), in the Rest Frame of the Front (RFF), i.e. where (fa^, = (1, 0, 0, 0), 
then P is unity. Otherwise, in the Rest Frame of the Gas (RFG), where = 
(1, 0, 0, 0), the escape probability P is -P*sc-ps = P^d^^^n Q{p^d?a^)/p^. 




Figure 2. A simple FO-surface in the Rest Frame of the Gas (RFG), where u'^ = (1, 0, 0, 0), 
including time-like and space-like parts. Then the normal vector of the FO front, d^cr^(a;), 
is time-like at the time-like part and it is changing smoothly into a space-like 4-vector in the 
space-like part. On these two parts of the hyper-surface, in the Local Rest Frame of the Front 
(RFF), d^aij,(x) points into the direction of the t' {x') axis respectively. 

The nominator depends on d^a^j,. In the following, we will take different 
typical values for d^cr^u, for characteristic regions (A, B, C, D, E, F) of the 
freeze-out hyper-surface (see Figure 2): 

■ A d^a^ = (1,0, 0, 0), leads to p^^d^a^ = P° > 0, as we have seen above 
for FO along the time axis in RFG, 

■ B = 7(a,6,0,0), 

■ C d^cTju = 7(1 + e, 1 — e, 0, 0), just above the light-cone, 

■ D d^a^ = 7(1 — e, 1 + e, 0, 0), just below the Ught-cone, 
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E 



7(c,(i,0,0), 



■ F Sa^ = (0, 1, 0, 0), leads to p'^d^a^ = p^, for FO along the spatial 
axis X in RFG, 

where the 7 factors serve the unit normalization of the normal vectors (see Figure 
2). The effect of this relativistically invariant FO factor leads to a smoothly 
changing behaviour as the direction of the normal vector changes in RFG (see 
Figure 2). 
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Figure 3. The contour plots of the escape probability factor as function of the particle momen- 
tum , -P(p), at different points of the FO hyper-surface are presented in six subplots. All plots are 
in the Rest Frame of the Gas (_RFG). ForregionA: d^a^ — (1, 0, 0, 0), and the probability fac- 
tor is one uniformly, for B: d^cr^ =7a(l,0.5,0,0),C(tirstrow)_D, d^cr^ =7,^(1,2,0,0) 
F: ct^a^ = (0, 1, 0, 0) (second row). The momenta are in units of [mc]. 

To calculate the parameters of the normal vector d'^a^ for different cases 
listed above, we simply make use of the Lorentz transformation. The normal 
vector of the time-like part of the freeze-out hypersurface may be defined as the 
local t'-axis, while the normal vector for the space-like part may be defined as the 
local x'-axis. As the d'^cj^ normal vector is normalized to unity its components 
may be interpreted in terms of and v^, as d^a^ = 7^(1, 0, 0), where 
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7ct = , ^ for time-like normals and 7cr = ^= — for space-like normals. 
This will lead to, 

■ A d?a^, = (1, 0, 0, 0), leads to P = 1, 

- B = 7^(1, va,0, 0), leads to P = 1^^^^^^, and = ^i=, 

■ C d^a^ = ^^(1 + 1 - 0, 0), leads to P = ^liiv^+v^+^Av'+v^) ^ 

m D dV^ = 7^(1-6,1 + 6,0,0), 

leads to P = 47g(P°+P-K^7..(P°-p-) ^ e(pMd3^^)^ 

■ E d^cr^ = 7(^(1,^0-, 0,0), 

leads to P = rf+HzE!) x GCp^^d^^^)^ 7. = ^7^, 

■ F d^o-^ = (0, 1, 0, 0), leads to P = |^ x eipi'd^a^), (see Figure 2). 

The resulting Phase Space escape probabilities are shown in Figure 3 for the 
six cases described above. 

In ref s. [10, 11 ] the post FO distribution was evaluated for space-like gradual 
FO in a kinetic model. Initially we had an equilibrated, interacting, PS distri- 
bution, fint{p, x), and an escape probability, similar to eq. (11), but simplified 
so that it depended on the angle of the two vectors only. After some small 
fraction of particles were frozen out as the FO process progressed in the front, 
the interacting component were re-equilibrated, with smaller particle number, 
smaller energy and momentum, to account for the quantities carried away by 
the frozen out particles. This was then repeated many times in small steps along 
the FO front and the frozen out particles were accumulated in the post freeze 
out PS distribution, ffree- The resulting distribution was highly anisotropic and 
obviously non-equilibrated. The details of the post FO distribution depend on 
the details of the escape probability, and on the level of re-equilibration of the 
remaining, interacting component. 

Bugaev assumed earlier [14], that the post FO distribution is a (sharply) "Cut- 
Juttner" distribution, but the above mentioned model shows that this can only 
be obtained if re-equilibration is not taking place. The kinetic model provided 
an asymmetric but smooth PS distribution [10], while the escape probability 
(11) yields a somewhat different, but also smooth PS distribution. These can 
be well approximated by the "Cancelling Juttner" distribution [15], Figure 4. 

Up to now we analysed the FO probability in the momentum space only. At 
this point we should also mention, that if we discuss gradual FO across a finite 
layer, the distance from a hypothetical or theoretical FO surface must also be 
part of the FO probability, because the probability of a particle has no more 
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Figure4. A simple FO-surface in the Rest Frame of the Gas (RFG) where ci^a,i = (1,0,0,0), 
including timelike and spacelike parts. Then d^a^, is t' = d'^a-^ for the timelike part, changing 
smoothly into x' = d^a-ft for spacelike normal vectors. 



collision on its way out depends on the local mean free path (m.f.p.) and on 
the distance from this surface. The m.f.p. is usually momentum dependent 
and this yields the secondary p— dependence of Pesc-ST{x;p). Of course the 
position of such a hypothetical FO surface has to be determined self consistently 
or within a realistic kinetic model. However, if we assume simultaneous FO 
and hadronization [16, 17] from supercooled plasma, there are no realistic and 
reliable dynamical models up to know, which we could use, so, a now the first 
self consistent approximate method is our only choice. 

4.2 Emission on an idealized hyper-surface 

If the m.f.p. or the coll. time are negligibly small compared to the ST 

dimensions of the system we can idealize the source distribution further as the 
source is shrinking to a hyper-surface (HS) in the ST, i.e.: 

S{t — Tpo) dr : time — like fronts, 



5{s — spo) ds : space — like fronts, 



dT 



27r Tcoll. 



exp 



(r - TFo) 



'^^coll. 
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and the integration in the normal direction to the front can be executed. Then 



Here e denotes integration over a 3 -dimensional HS. If all particles described 
by fe.q. are emitted during the process, Ci = 1. If we do not perform the 
momentum integral we recover the modified "Cooper-Frye" freeze out formula 



Notice that we have changed the notation of the PS distribution from f^q to fpo 
denoting the post FO distribution integrated across the front. 

Eq. (13) simplifies the process to a great extent, especially if we consider 
that feq changes during the detonation deflagration or freeze out process as 
we go across the front self-consistently. In a dynamical process across the 
front, even if /eg is always an equihbrated distribution, its parameters change 
in a wide range [10, 11]. In the FO process the post FO distribution, after the 
process is completed is very far from any statistically equilibrated distribution. 
This should be taken into account if we apply the assumption of an idealized 
FO surface. The determination of the FO surface is an involved task, both 
physically and technically. In complex modeling tasks we perform the self- 
consistent determination of the surface in simplified models. In these models 
we can pinpoint the FO conditions and evaluate the post FO distributions. 

Then based on this knowledge we determine the FO surface in the large 
scale, 3+1 dimensional system by analysing the complete ST history of the 
fluid dynamical stage calculated. This demanding task was completed recently 
by BemdR. Schlei (of Los Alamos) [20], and this enables us to visualize the 
ST evolution of the surface where the emitted particles originate from. 

5. Two particle correlations 

In case of two particle correlations we have to return to eqs. (8) or (9), because 
the emission point in ST is crucial. Although, eqs. (12) or (13) also contain a 
ST dependence, this is reflecting the properties of the idealized hypersurface, 
where the process is taking place. This approximation is justified only if 

■ The ST domain of the emission is a negligibly narrow layer and we are 
sure that the internal processes within the front are not relevant. 

■ If we use a realistic estimate for the post FO distribution in these equa- 




(12) 



[18, 14, 19] 
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The PS density of created (or frozen out) particles, Ncr, is given in terms of 
a "Source function", S{x,p), as [21] ^ 

^ = 1 Six,p) d'x , (14) 

which can be compared to eqs. (8,1 1), so that 

j d'^X feqip,x) Pesc-PHip;x) Pesc-ST{x\p) Or (15) 



dNcr 
(Pp 



, . 

d'^X feq{p,x) Pesc-PH{p\x)/p^ Pesc-ST{x\p) . (16) 



S(j>;x) 

Thus, assuming emission or freeze-out from an interacting gas the source 
function is 

S{x,p) = feq{p, x) Pesc-PH{P\ x) Pesc-ST{x;p) , (17) 

here we assumed that the emission or escape probability is factored to a PS and 
a ST dependent part. 

Two particle correlations provide information about the location of FO points 
in the ST, providing the possibility of deeper insight into the reaction mech- 
anism. This obviously means that there is an essential difference between 
reaction models assuming (i) FO in an extended ST layer or (ii) FO accross an 
idealized ST hypersurface. The latter assumption is obviously not realistic. 

Nevertheless, the idealized hypersurface assumption may not be so bad as 
one would think. Even in this case parts of this surface can represent early 
times while others late emissions. The spatial locations of different parts of 
the surface can be far from each other. If the ST dimensions of the idealized 
FO hypersurface are much larger than the thickness of the real FO layer, the 
hypersurface assumption may be satisfactory. In any case, two particle cor- 
relations are much more sensitive to this simphfying assumption than single 
particle data. 

Note also that the thickness of the FO layer is dependent of the hadronic 
species measured. Thus, the hypersurface idealization is applicable for high 
multiplicity hadrons, protons and pions, while particles with low multiplic- 
ity, low cross-sections or low creation rates are not described well with this 
approximation. 

6. Simultaneous Hadronization and Freeze-Out 



If one assumes hadronization in thermal and fluid dynamical equilibrium 
via homogeneous nucleation [22] or similar processes, this leads to a lengthy. 
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gradual hadronization and freeze-out. Such a long-lived particle source should 
be detected by pion interferometry, but recent data from RHIC ruled out this 
scenario. Therefore we do not assume a lengthy thermalization and chemical 
equilibration process at freeze-out either. As we described in ref. [23], we 
assume simultaneous FO and hadronization at the end of the heavy ion reaction. 

We assume that pre-hadrons or quasi-hadrons are formed in the cooling 
and expanding plasma before we reach the critical temperature and density, 
and especially when the plasma supercooles. Supercooling is possible even 
if the phase transition is just a smooth, but sharp cross-over, which is always 
the case in small, finite systems. Then hadrons are formed out of thermal 
equilibrium, which are frozen out (i.e. never collide) after their formation. No 
detailed hadronization models exist for this kind of mechanism. The closest 
are coalescence or recombination models, which usually describe the hadron 
abundances well. 

As discussed in [23] and [24], the similarity of the results arising from sta- 
tistical (thermal) and coalescence models, is due to the fact that the Canonical 
Ensemble (CE) is generated if the average energy of the system is the same in 
the elements of the ensemble and the states included are generated with equal 
a priory probability. For most hadrons these conditions are satisfied, so the 
two models yield similar results and therefore we use the statistical model to 
estimate local hadron abundances at the end of the reaction [23]. 

The differences arise for hadrons where the assumption of equal a priory 
formation probabihty does not hold. An example is the A(1520) resonance, an 
excited state where the radial wave function includes both s- and p- waves, and 
so the formation cross-section is much less than for other hadrons. Therefore 
the abundance of these hadrons is much smaller, and this is well described by 
coalescence type models, while statistical models fail to reproduce the smaller 
abundance (as the cross-section does not appear in these models). 

For our purposes, however, the statistical model is adequate, because it in- 
cludes all necessary statistical weights, (so it works well for the high multi- 
plicity hadrons,) while low multipUcity resonances are not reproduced by fluid 
dynamical models anyway. 

7. Summary 

Recent experimental data indicate that detailed fluid dynamical data are be- 
coming available at RHIC energies. This will provide us with the possibility to 
test the QGP Equation of State in heavy ion reactions. This will require further 
experimental efforts that provide us with the complete spectrum of collective 
flow data. To draw the right conclusion will also require a detailed and realistic 
theoretical reaction model which can simultaneously describe single particle 
and two-particle observables. 
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Notes 

1. Note that in some cases the source function is defined differently [25], so that the p*^ factor is not 
included in the denominator: = / S'{x,p)d*x . 
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